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Several models of neutrino masses predict the existence of neutral heavy leptons. Here, we re¬ 
view current constraints on heavy neutrinos and apply a new formalism separating new physics 
from Standard Model. We discuss also the indirect effect of extra heavy neutrinos in oscillation 
experiments. 


I. INTRODUCTION 

In the Standard Model (SM), neutrinos are massless 
particles contradicting the experimental observation of 
neutrino oscillations, hence physics beyond the SM is re¬ 
quired. 

An effective lepton number violation dimension- 
five operator O 5 oc LL$< 1 >, can be added to the 
SU(3)c G SU(2)l 0 U(1)y model in order to introduce 
neutrino masses mm, where L is one of the three lepton 
doublets and is the SM scalar doublet. 

After electroweak symmetry breaking, Majorana neu¬ 
trino masses are induced, being proportional to ($)^ and 



FIG. I: Dimension five operator responsible for 
neutrino mass. 


implying lepton number violation. Hence, the smallness 
of neutrino mass, compared to the masses of the SM 
charged fermions, arises from the smallness of the co¬ 
efficient in front of the operator O 5 associated with the 
lepton number violation by two units (AL = 2). 

Unfortunately we cannot say too much more about this 
operator. We do not have any clue about its mechanism, 
nor its mass scale, nor its flavour structure. 

A common possibility is to assume that O 5 is induced, 
at the tree level, by the exchange of heavy messenger 
particles. In this way, seesaw models postulate neutral 
heavy states act as messenger particles to induce neu¬ 
trino mass. For instance, a right-handed neutrino could 
be included associated to each of the three isodoublet 
neutrinos (Type I seesaw) [3]: 


rrii, = Ao 


M 


( 1 ) 


The existence of processes with AL = 2, such as neu¬ 
trinoless double beta decay, or lepton flavour violation 
processes (LFV) as ^ ey , would give hints on the pos¬ 
sible existence of these heavy Majorana neutrino messen¬ 
gers. Hence we could find signatures of heavy neutrinos 
and their mixings by studying this kind of processes. 


II. THE METHOD 
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As we said in the previous section, heavy neutrinos are 
introduced in several extensions of the SM such as linear 
and inverse seesaw models mm, leading to a rich struc¬ 
ture in the lepton mixing matrix. In order to work with 
this kind of models we will use a symmetric parameteri¬ 
zation, consistent with the general formalism [4], neatly 
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separating “new physics (NP)” and “Standard Model 
physics (SM)”. 

For the case of three light neutrinos and n — 3 extra 
heavy states, we can construct the mixing matrix U as the 
product of LUij rotation matrices (Okubo’s notation [U 
0): 

jrnxm _ , , . , , , 

U — UJn—ln^n—2n n 

Wn-2n-l . . . 1^2 3 3 Wi 2 ( 2 ) 


In summary, by choosing a convenient order for the 
products of the rotation matrices, we can obtain a 
parameterization which puts all the information in a con¬ 
venient form. 

A. Simplest extension of SM: 3 -|-1 neutrinos 

The formalism for the simplest extension of the SM 
includes one extra right handed singlet 


/ Cij 0 e 

Wy = 0 1 0 ( 3 ) 

0 dj j 

where = cos 9 ij and Sij = sin 6 ij. 

The mixing matrix U can be decomposed in a new 
physics part and its Standard Model part 


jj-nxn 

— R^^ R^^ 


( 4 ) 

rSM 

~ ^n—ln^n—2n • 

= 012 3 ^13 ^12 

■ ■ n ■ 

■ ■ a;3 4 012 4 a;i 4 ( 5 ) 

(6) 


and it can be divided in four blocks 


U = 



( 7 ) 


where N is the block corresponding to the standard three 
neutrinos, including their mixings between them and the 
extra neutrinos. 

At the same time, the matrix N can be also decom¬ 
posed as 


with the mixing relations between the gauge and mass 
eigenstates given as m 

3 

—'^^Wka^aL + (H) 


The unitary mixing matrix {7^x4 written as 


jjAxA _ 


' j\^ 3 x 3 5 * 3 x 1 ' 

2 ^ 1 x 3 ^1x1 



Ne2 

A^e3 

5e4\ 

N^i 



Sfx A 

Nrl 

Nr2 

A'r3 

SrA 

yT 4 i 

T42 

743 



( 12 ) 

where is again the sub-matrix related with the 

standard neutrinos Eq.Q It is important to notice that 
jg unitary whereas is unitary because in¬ 

cludes all neutrinos in the model. 

Comparing the terms in with the terms of 

we obtain the following expressions for the a factors 


N = 


0 0 \ 
021 022 0 1 
\03i 032 033/ 


( 8 ) 


where is the usual mixing matrix of the Standard 

Model and the matrix includes all the new physics 

information through the o^- parameters (a more com¬ 
plete discussion is given in m) 


Oil = Ci 4 , 

O22 = C24 , 

O33 = C34 , 

O21 = S24 Si4 , 

O32 = S34 S24 , 

O3I = S34C24S14. 


( 13 ) 


On 

— Cl n 

C In— 1 Cl 

n —2 ■ 

• ■ Ci4 


O 22 

~ C 2 n 

C2n-lC2 

n-2 ■ 

• ■ C 24 


033 

= C3n 

C3n-lC3 

n-2 ■ 

• ■ C 34 

(9) 

021 

— C 2 n 

C2n-1■■ 

• C 25 ■ 

S 24 S 14 



+ C2n 

. . . C26 S 25 S 15 

Cl4 “t” ^2nSln Cln—1 Cin—2 ■ 

■ ■ Ci4 


B. Aplication to 3 +3 model 

Usually, more than one extra neutrino is introduced 
in the theory, as in sequential-type seesaw mechanisms 
where 3 (Type I) or 6 (Inverse and Linear) extra singlets 
are included with the SU{ 2 )l SM doublets. 
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For such models our parameterization becomes 


jj6x6 


( jY3x3 g3x3 
2^3x3 ^3x3 


(14) 


with these expressions for a the parameters 


ail = Ci6 Ci5 Ci4 , 

022 = C26 C25 C24 , 

033 = C36 C35 C34 , 

O2I = S26S16C15C14 + C26S25S15C14 

+ C26 C 25 S 24 Sl4 , (15) 

O32 = C36 C35 S34 S24 + C36 C35 S25 C24 
+ S36 S26 C25 C24 , 

O3I = C36 C35 C34 S34 C24 Si4 + C36 S35 C24 Sl5 Cl4 
+ S36 C26 S16 C15 Ci4 + C36 S35 S25 S24 Sl4 

+ S36 S26 C25 §24 §14 + §36 §26 §25 §15 C14 ■ 


For the electron (anti) neutrino survival probability, we 
get the expression 
3 

Pee = ^ |lVe,f |iVe,f 

3 

k>j 


4P 


(18) 


and using Eq.pd|) and Eq.p^ we obtain 
4 


Pee = O? 


1 — cos^ 013 sin^ 2012 sin^ 


/ Am^yL 
( 4P 


— sin^ 2013 sin' 


Am\^L 

4E 


(19) 


Note that the effect of the extra neutrinos is totally in¬ 
cluded in the afi factor, illustrating the utility of our 
symmetric formalism. 

Considering now only one extra fourth neutrino, 


Eq.(18) would change to 

cos'‘0i4 = (1 - l'S'e4|^) 


Pe, 


( 20 ) 


III. OSCILLATION CONSTRAINTS 

The general expression for the survival and conversion 
neutrino probability is given by m 

Pc.p = (5„;3 - 4 ^ me Upk Pai Sin^ 

k>j 

( Am? L\ 

( 16 ) 

/ 

where Sap appears due to the unitarity of the mixing 
matrix. However, in a model with extra heavy neutri¬ 
nos, this equation will change because the mixing ma¬ 
trix describing the three standard neutrinos (A^^^ in our 
symmetric notation) will not be unitary and the effective 
probability will not be normalized to 1. In this case, 
the probability includes the W terms from the truncated 
matrix 

3 

Pcp = Y.K.NpkN^.N^p^ (17) 

j,k 

- AY,^z[Nl^NpuN^,Nl] sin- 

k>j 

+ 2 ^ am Npk N^j N*p^] sin 

k>j 


One could be tempted to study whether the presence 
of an extra light singlet leptons (sterile neutrinos) could 
play some role in the reported neutrino anomalies (as 
MiniBooNE m)- Unfortunately this is not the case and 
we stick to the (natural) assumption that the extra states 
are heavy and do not take part in oscillation effects. In 
this case one can get a constraint from some reported 
combined analysis |13j 

sin^ 014 = 15^41^ < 0.04 (90%C.L.) (21) 

ail = cos^ ^*14 < 0.96 (22) 

IV. FUTURE OSCILLATION EXPERIMENTS 

As we can see from the previous section, it is not pos¬ 
sible with current experiments to obtain a very strong 
constraint on the new physics a parameters so we con¬ 
sider future experimental proposals such as LENA El- 

LENA is a future neutrino experiment which will use a 
^^Cr artificial neutrino source with 5MCi intensity, pro¬ 
ducing a total of 1.9 X 10® neutrino events. The expected 
number of neutrino events for an energy recoil of the elec¬ 
tron in the range from 200 to 550 keV in the presence of 
an extra heavy neutrino would be given by 

A, = af 1 n, c^cr At J (T, T') dT dT' (23) 



'Aml^L" 

AE 


Arr^ 

2E 
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FIG. 2: Estimation of LENA sensitivity. 


where rie is the number of electron targets, (j)cr is the 
neutrino flux coming from the source, At is 28 days which 
corresponds to the half-life of the source, and R {T, T') 
is the resolution function 


i?(T, T') 



exp 


[T-T'f 

2 ct 2 


(24) 


where T is the recoil energy, T' is the true energy and 
cr = 0.075-\/r/MeE is the expected energy resolution. 

An estimate of the expected sensitivity as a function 
of the total percent error of the experiment can be per¬ 
formed in advance, being shown in the Fig. 


V. OTHER CONSTRAINTS 


also the excluded region from neutrinoless double beta 
decay experiments |24j , valid only if the heavy neutrinos 
are Majorana particles. For completeness we show the 
bounds on the mixing with muon and tau neutrinos in 
Fig.g 

From a more complete combined analysis of these data 
we can get constraints on our new physics parameters 
associated to the presence of the heavy neutrinos 


VI. CONCLUSIONS 

Extra neutral heavy leptons are motivated in order to 
introduce neutrino mass but no positive evidence of these 
particles has been found so far. 

Signatures of these heavy neutrinos arising from their 
mixings with the light ones could be searched at labora¬ 
tory experiments. The study of these bounds would be 
useful to shed light upon the mass generation mechanism 
of neutrinos and probe the scale of new physics. 

Models beyond SM, such as seesaw models, imply a 
very large number of parameters. The symmetric param- 
eterizatio of the neutrino mixing matrix, describing the 
charged current, provides a very useful way to separate 
new physics from SM effects, concentrating the informa¬ 
tion and making easier to work with it. A more detailed 
account of our work will be described elsewhere pS] . 


The effects of heavy neutrinos would show up as peaks 
in the leptonic decays of pions and kaons or from their 
direct production at higher collider energies [S] . One can 
perform an analysis of all these data, and combine the 
corresponding restrictions on heavy neutrinos parame¬ 
ters. 

Fig. [^compiles the bounds on the heavy neutrino mix¬ 
ing with electron neutrinos at 90%C.L. It agglutinates 
bounds coming from peak searches at lepton decays, as 
TT —)■ e i/ [m [IS] and K ^ ev m HH); meson decays 
at PS191 [HI, NA3 [20| and CHARM [H] and de¬ 
cays at DELPHI [22| and L3 [23] . In Fig. we show 
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